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Abstract. We study the sensitivity of the sound speed 
to relativistic corrections of the equation of state (EOS) 
in the standard solar model by means of a helioseismic 
forward analysis. We use the latest GOLF/SOHO data 
for £ = 0,1,2,3 modes to confirm that the inclusion of 
the relativistic corrections to the adiabatic exponent T\ 
computed from both OPAL and MHD EOS leads to a 
more reliable theoretical modelling of the innermost layers 
of the Sun. 
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1. Introduction 

It has recently been shown (Elliott & Kosovichev 1998) 
that the inclusion of relativistic effects in the equation of 
state (EOS) leads to a very good agreement between the 
solar models and the seismic Sun. In particular, the inver- 
sions of SOI-MDI/SOHO p-mode frequencies for the adi- 
abatic exponent Ti show that MHD EOS reproduces the 
interior of the Sun with great accuracy, when the relativis- 
tic contribution to the Fermi-Dirac statistics is included. 
It is thus interesting to approach the same problem 
by means of the forward analysis by comparing the the- 
oretical eigenfrequencies with the observed ones. Unfor- 
tunately this method is not directly applicable since our 
description of the outer layers of the Sun is still far from 
complete and many theoretical uncertainties would influ- 
ence our conclusions. However, since such small effects in 
solar EOS are most important only in the deep interior, it 
is possible to make use of the acustic mode frequency small 
separation diagnostic, 5vi^ n = Vi yTl — z^+2,n-i, f° r spherical 
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harmonic degrees £ = 0, 1 and radial order n 3> £ (Tas- 
soul 1980). The main property of this quantity is that it 
is strongly sensitive to the sound speed gradient near the 
solar centre while it is weakly dependent on the details 
of the treatment of the outer layers. Since the relativistic 
effects manifest themselves mainly through a depletion of 
0.1% — 0.2% of the adiabatic index Ti, we expect a quan- 
titatively similar change of sound speed gradient in the 
solar core. 

The acoustic mode frequency small separation analysis 
has recently been used for estimating the seismic age of 
the Sun (Dziembowski et al. 1998) and the related impli- 
cations of the uncertainties in the Sn astrophysical factor 
determinations (Bonanno & Paterno 2001). 

Here we show that the mentioned above analysis can 
also be used to verify how the different physical character- 
istics of the MHD and OPAL EOS reflect on the accuracy 
of the description of the stratification of the internal lay- 
ers of the Sun. On performing a x 2 analysis of the latest 
published GOLF/SOHO data for different solar models, 
we confirm the main conclusion of Elliott & Kosovichev 
(1998), based on an inversion analysis, that the inclusion 
of the relativistic effects in the EOS is in any case required 
to improve the accuracy of solar models, independent of 
which EOS is used. 

2. The solar model 

In our analysis we used the GARching SOlar Model 
(GARSOM) code which has been described in detail in 
Schlattl et al. (1997). It includes the latest OPAL-opacities 
and either OPAL or MHD EOS, and it takes into account 
the microscopic diffusion of the elements heavier than hy- 
drogen. Our standard solar model has been verified in de- 
tail in Turck-Chieze et al. (1998) and found in good agree- 
ment with other up-to-date solar models, and, in particu- 
lar, it is consistent with the observed L Q and T e ff within 
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10~ 4 , at an age of 4.60 Gy, adopting the surface value 
Z/X = 0.0245. We then included the relativistic correc- 
tion leading term to the adiabatic index Ti derived from 
the relativistic evaluation of the Fermi-Dirac integrals of 
the EOS in the solar core by means of the expression (El- 
liott & Kosovichev 1998): 
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where T is a dimcnsionless temperature in units of m e c 2 /fc, 
with m e the electron mass, c the light speed in vacuum, k 
the Boltzmann constant, and X the hydrogen abundance 
by mass. As expected, the correction to Ti is negative, 
namely IY ro i < ri. nr , since Ti tends to shift from the 
non-relativistic value of 5/3 to the extremely relativistic 
one of 4/3. 

The corresponding relativistic corrections to the lead- 
ing terms for sound speed, c s , and density, g, are respec- 
tively: 
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where tp is the degeneracy parameter, that is about -1.14 
at the Sun's centre, and decreases noticeably toward the 
surface, the partial degeneracy being completely removed 
at 0.4 R . 

The behaviour, as functions of the fractional radius, of 
the relative differences between the quantities Ti, c s and g 
calculated with relativistic corrections and without them 
is shown in Fig. 1. 

The term -(15/64\/2)f e^ = <5P/P - Sg/g, in Eq.(2) 
is negligible with respect to 5Ti/T\ indicating that the 
relativistic corrections to the pressure, P, and density, g, 
cancel each other almost completely and the correction to 
c s is entirely dominated by the correction to T\. Also the 
term e^/4\/2 in Eq.(3) is negligible with respect to unity, 
indicating that in the solar case the coupling between de- 
generacy and relativistic effects is weak. 




Fig. 1. Behaviour, as functions of the fractional radius, of 
the relative differences between relativistic and non-relativistic 
quantities <5y/y = (y r8 i — ynr)/y nr , where the ys stand for Ti 
(continuous line), c s (dashed line), and g (dashed-dotted line) 
respectively. 



where i stands for £ = 0, 1, m = 10 and M = 26. Fig. 2 
and Fig. 3 show the behaviour of the terms in the sum 
defined in Eq.(4) in non-relativistic and relativistic cases 
for MHD and OPAL EOS respectively. The difference 
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3. Results with GOLF/SOHO data 

We used the latest GOLF/SOHO data for I = 0,1,2,3 
obtained with long time series and by taking into account 
the asymmetric line profile in data reduction (Thiery et 
al. 2000). In particular, we determined the acoustic mode 
small spacing difference bv^^ for I = 0, 1 and n>( for our 
solar model, and studied the difference <5^, n ,0 — <5^i. n , model 
between data and model. We then constructed the two x 2 
indicators (Dzicmbowski et al. 1998, Schlattl et al. 1997) 
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Fig. 2. Relativistic (continuous line) and non-relativistic 
(dashed line) contribution to the x 2 calculation for MHD EOS. 



between the relativistic and non relativistic case is larger 
for £ = in the frequency range beetwen 2000 and 2500 
//Hz, and for £ = 1 beetwen 2500 and 3000 //Hz. The \ 2 
results are shown in Table 1 where it is possible to note 
that the models with relativistic corrections have rather 
smaller x 2 s and there is no significant difference between 
Xq and x 2 calculated for OPAL and MHD EOS. However, 
MHD EOS appears to be slightly favoured with respect to 
OPAL EOS. 
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Fig. 3. Rclativistic (continuous line) and non-relativistic 
(dashed line) contribution to the \ 2 calculation for OPAL EOS. 



Table 1. x 2 results in non relativistic (NR) and relativistic 
(REL) cases for MHD and OPAL EOS. 



EOS xd(NR) Xo(REL) 



*f(NR) Xi(REL) 



MHD 
OPAL 



1.73 
1.91 



1.41 
1.41 



2.13 
2.32 



1.67 
1.94 



4. Conclusions 

Our results show that the acoustic mode frequency small 
separations are sensitive to the inclusion of the rclativistic 
effects. It would be interesting to discuss the relevance of 
these effects in the helioscismic determination of the solar 
age and related problems with Sn uncertainties. We plan 
to address this issue in a forthcoming communication. 
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